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Abstract — A linear output feedback control 
scheme is developed for a coupled map lattice system. 
Hoc control theory is used to make the scheme local: 
both the collection of information and the feedback are 
implemented through an array of locally coupled con- 
trol sites. Robustness properties of the control scheme 
are discussed. 

I. Introduction 

Learning to tame spatiotemporal chaos in spatially ex- 
tended nonlinear systems is very attractive due to a 
large number of potential applications. Some of these 
are continuous, such as turbulence jy, plasma insta- 
bilities H and chemical reaction systems || , some are 
discrete: neural networks and distributed memory 
systems are only a few examples. The main objective 
is usually to stabilize some suitable unstable periodic 
orbit (UPO), or a group of orbits, embedded in the 
chaotic attractor of the system. 

Although spatially extended homogeneous systems 
could be treated as a special case of the high- 
dimensional chaotic systems, some of the practical is- 
sues, that arise in the control problem are quite specific 
and could be best handled by taking into account the 
spatiotemporal structure of the system and the con- 
trolled state in general and their symmetry properties 
in particular ||. 

In the present paper we will illustrate the control al- 
gorithm applying it to the general coupled map lattice 
(CML), originally introduced by Kaneko j6|: 
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/(4) + e(/(4-i)-2/(4) + /(4+i)), (i) 



and considered to be one of the simplest models, pos- 
sessing the essential properties of an extended spa- 
tiotemporally chaotic system. 

There are many ways to achieve the stabilization of 
a non-chaotic trajectory. However, the requirements 
imposed by different control algorithms and their per- 
formance could vary widely. For instance, it was shown 
@, that a number of UPOs of the CML (|l|) could be 
stabilized with feedback applied through a periodic ar- 
ray of controllers. Although limited knowledge of the 
system state was required, the density of controllers 



had to be extremely high for the control to work. Re- 
arranging the controllers, one can significantly reduce 
their density and improve the robustness characteris- 
tics of the control scheme § at the expense of requir- 
ing additional information about the system state. In 
the present paper we will show how the CML can be 
controlled using low density of controllers and requir- 
ing very limited information about the system state. 

II. The system 

Rewrite eq. (jl]) , adding to it the uncorrelated random 
noise (wjwj, ) = a 2 5tt'Su' and applying control pertur- 
bations ut — Gfc(z', z* _1 , • • •) at sites ik, k = 1, • • • ,m: 

4 +1 = eftt-i) + (1 - 2e)/(4) + tf(4fi) 
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assuming, that the lattice is finite, i = 1, 2, • • • , n, and 
periodic boundary conditions z\ +n — z\ are imposed. 

Due to the translational symmetry of the CML ([!]) 
additional parameters can only enter the evolution 
equation through the nonlinear local map function, 
which we choose as f(z) = az(l — z), emphasizing, 
that the only result affected by this particular choice 
is the set of existing periodic trajectories. In partic- 
ular, for any choice of f(z), the homogeneity of the 
system response to the perturbation of any internal 
parameter (a and e in our case) makes it impossible to 
use either internal parameter for control. 

Linearizing equation (Q) around the period-r tar- 
get UPO z , z 2 , • • • , z T and denoting the displacement 
x* = z — z , we obtain 
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(3) 



where A\, — djZ* + (z*) is the Jacobian and the ma- 



trices B l N 



and B] 



l3 - Efc specify the 

response of the system to the external noise w\ and 
the applied feedback u\ (also called the input). 

Finally, assume that only q functions r\\ = i?i(z*) 
(called the output) of the system state are accessible 
to measurement. Denoting CL = djHiffi) we obtain 
for the linearized output: 



(4) 



III. The control scheme 

The algorithm presented below allows one to deter- 
mine whether the feedback u* stabilizing the chosen 
UPO can be obtained as a function of the output y*, 
and determines the solution, which minimizes the noise 
amplification factor or induced-power-norm 
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where the r-dimensional performance vector 



a* = C 



(5) 



(6) 



gives the deviation of the system from the target state, 
and the power norm is defined as 
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(7) 



The solution to the time-periodic output feedback 
problem (d^H) c an be obtained using the generaliza- 
tion of the results of Hoc control theory || for linear 
time invariant (LTI) systems. In particular, Dullerud 
and Lall have shown fTo|j , that if a locally stabilizing 



linear feedback u exists, it could be written as 
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= Cr 



(8) 



where Aq, B l Cl C l c and D l c are matrices with the same 
periodicity r as the target orbit z*, and v* is the p- 
dimensional internal state of the controllers. The stan- 
dard state feedback law u* = K 'x* used in Q , is seen 
to be just a special case of this general setup. 

Construct constant block diagonal matrices A,B,C, 
Bn,Cn and Dn according to the following rule: 
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(9) 
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For r > 1 define a rn x rn cyclic shift matrix 
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In the time-invariant case (r = 1) set Z = I n xn- Also 
introduce the notations Q > for positive definite, 
Q > for semi-positive definite matrices and Q f° r 
the transpose of Q. 

It can be shown |l(J, that a stabilizing solution (||) 
with p > n such that 7 < 1 for the system (||-||) exists, 
if and only if there exist block-diagonal matrices R > 
and S > 0, satisfying 
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where Pr, P5, Or and Os are given by 
'A^Z^SZA-S A^Z^SZB N 



Ps = 



Pr = 
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(12) 



and the unitary matrices Nr and Ns satisfy 
Im Nr = ker [ fit 



ImjVs=ker[C* pXn ] . 



To minimize 7, rescale C l N and D l N , 



(13) 

such that the 
above condition tests for 7 < 70 instead of 7 < 1 and 
decrease 70 until the test fails; standard software exists 
to do this. If there is any linear stabilizing controller, 
we can therefore find it using this algorithm. 

If R = diag(i? 1 , . . . , R T ) and S = diag(S' 1 , . . . , S T ) 
are determined, one can find the matrices in (|J) using 
the following procedure. First, construct nonsingular 
matrices M t and N t , such that 

M t Nj = I-R t S t . (14) 

Determine the matrix X t as the unique solution of 



Next, define the matrices 
A 1 nx 



A t 
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and then define 
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Finally, the matrices A^^B^jC^ 
tracted from the solution 

Jt 



and D l c 



A 1 

r l 



to the linear matrix inequality 

H t + Q\j}P t + P}JtQt <0. 



(17) 
are ex- 

(18) 
(19) 



Linear matrix inequalities (LMI) like ( JlT| ) and ( ^9[ ) 
can be conveniently solved using the tools of convex 
optimization theory. The big practical advantage of 
this technique is the guaranteed convergence. 
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Figure 1: Stabilization of unstable periodic orbits of 
the noisy system: (a) S1T1, (b) S8T2 and (c) S8T4. 
4th lattice variable z\ is plotted. The time of capture 
is taken to be t = 0. The noise strength is a = 10~ 6 . 



IV. Control of large lattices 

Although, using the above algorithm, we can in prin- 
ciple obtain the stabilizing feedback (|^) for a system 
(^|) of arbitrary size, solving matrix inequalities involv- 
ing large matrices requires considerable computational 
resources. 

This problem could be avoided using distributed 
control approach. The idea is to subdivide the com- 
plete system into a number of weakly interacting sub- 
systems, and learn to control each of the subsystems 
independently, neglecting interactions with other sub- 
systems. Finally, the control can be adjusted to cor- 
rect for interactions by introducing coupling between 
formerly independent controllers. 

Since the coupling in our model is local, we can par- 
tition the whole lattice into a number of identical sub- 
domains of length n p <C n, each interacting with two 
adjacent subdomains. The original problem is thus 
reduced to the problem of controlling an isolated sub- 
domain of limited length n p (we drop the index below). 
We impose periodic boundary conditions on each sub- 
domain to allow the existence of unstable orbits peri- 
odic in space as well as time. 

The symmetry properties of the CML (^) determine 
|5j , that the minimal number of controllers required is 
two. Placing them at the boundaries of the subdomain 
allows one to change the boundary condition at will, 
as well as correct for interactions between adjacent 
subdomains, by adding appropriate perturbations to 
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Figure 2: Stabilization of the homogeneous time- 
period-11 (S1T11) UPO. The time of capture is taken 
to be t = 0. The noise strength is a = 10 -6 . 



the feedback 



This defines the matrix 



-B* — SuSjl + Si n Sj2- 



(20) 



In order to calculate these perturbations we will 
have to introduce coupling between controllers of ad- 
jacent subdomains. Specifically, we will need to ex- 
change the information about the state of the system 
in the neighborhood of the boundaries (and therefore 
controllers), i.e. at least the variables x\ and x f n should 
be measurable. This defines the minimal realization of 

= 2, which we use below. 



the matrix C* 



B%,q 



V. Comparison of H2 and Hoo approaches 

In order to compare the results of the proposed ap- 
proach with those, obtained using linear quadratic 
(H2) theory for the state feedback ||, we select a sim- 
ilar optimization criterion. Specifically, we take 
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(21) 



z* = {x';u*} and, consequently, 



lim 
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(22) 



We demonstrate the approach by stabilizing a 
number of UPOs of the noisy CML (||) with n = 8 
sites, a = 4.0 and e = 0.33. The feedback (||) is cal- 
culated using the algorithm outlined above. Figure 
^ shows the process of capturing and controlling the 
steady homogeneous state (S1T1), the time-period-2 
space-period-8 (S8T2), and the time-period-4 space- 
period-8 (S8T4) orbits. 

The real power of the approach, however, can 
be full appreciated only in application to orbits of 
very high periodicity, where the accurate treatment 
of the effects of noise is of ultimate importance. Any 
method based on the reduction of periodic trajectories 
to steady states will fail for orbits of sufficiently long 
period. The approach does not suffer from this 
limitation. Indeed, we have observed stabilization of 
a number of periodic orbits with period t > 10. One 
such example is presented in Fig. 0. 




Figure 3: The largest length of the lattice, which 
can be controlled with two pinning sites: the triangles 
represent the length obtained using Hoc control, the 
circles show the data, obtained using H2 control in 
Ref. and the curves show theoretical estimates 

The noise strength is a = 10 -14 and 



23J) and (|25j) 
a = 4.0. 



Noise limits our ability to control arbitrarily large 
systems with local interactions, using just two con- 
trollers. Rather simple arguments show H, that the 
size of the largest system, that could be stabilized in 
the presence of random perturbations w , could be es- 
timated using the controllability condition, if complete 
information about the state of the system is available. 
So, for a steady uniform state one obtains 



e > 0.5 

2(log(e)-A roox )- 1 log(V), e<0.5, 



(23) 



where X max is the maximal Lyapunov exponent. 

If however only partial information about the state 
of the system is available, additional requirements ap- 
pear. Any control algorithm utilizing output feedback 
essentially consists of two major stages: observation 
and control. During the first stage information about 
the system is collected and processed to recreate the 
state of the system. During the second stage, control 
perturbations are applied to bring the system to the 
desired state. As a result the control scheme should be 
able to tolerate uncertainties introduced during both 
the control and the observation stage. 

The estimate (^3|) reflects the requirements imposed 
by the control stage. Additional requirements, intro- 
duced by the observation stage can be similarly esti- 
mated using the observability condition, 



rank[ct 



(24) 



which determines whether the state of the system can 
be extracted from the observed data ([!]), and for C = 
coincides with the controllability condition. 



Careful consideration shows, that the addition of the 
observation stage effectively doubles both the length of 
the control cycle and the length of the lattice. As a 
result, the maximal length of the system, that can be 
successfully stabilized using Hoe control is halved: 



,((t) =na(tr)/2. 



(25) 



The maximal length n 00 (cr) can be obtained numer- 
ically by choosing the fixed point as the initial con- 
dition and monitoring the evolution of the system in 
the presence of noise under control (||). The results 
are presented in Fig. ||. One can see that the estimate 
(25) approximates the actual results rather well. 
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through grant no. DMR-9013984. 

References 

[1] C. Lee, J. Kim, D. Bobcock and R. Goodman, "Ap- 
plications of neural networks to turbulence control 
for drag reduction," Physics of Fluids, V 9, N 6, 
pp.1740-1747, 1997. 

[2] A. Pentek, J. B. Kadtke and Z. Toroczkai, "Sta- 
bilizing chaotic vortex trajectories - An example of 
high-dimensional control," Phys. Lett. A, V 224, N 
1-2, pp.85-92, 1996. 

[3] V. Petrov, M. J. Crowley and K. Showalter, 
"Tracking unstable periodic orbits in the Belousov- 
Zhabotinsky reaction," Phys. Rev. Lett., V 72, N 
18, pp.2955-2958, 1994. 

[4] C. Lourenco and A. Babloyantz, "Control of spa- 
tiotemporal chaos in neuronal networks," Int. J. 
Neu. Sys, V 7, N 4, pp.507-517, 1996. 

[5] R. O. Grigoriev and M. C. Cross, "Controlling- 
physical systems with symmetries," submitted to 
Phys. Rev. E. 

[6] K. Kaneko, "Period-doubling of kink-antikink pat- 
terns, quaziperiodicity in antifcrro-like structures 
and spatial intcrmittcncy in coupled logistic lattice 
- towards a prelude of a field theory of chaos," Prog. 
Theor. Phys., V 72, N 3, pp.480-486, 1984. 

[7] G. Hu and Z. Qu, "Controlling spatiotemporal 
chaos in coupled map lattice systems", Phys. Rev. 
Lett., V 72, N 1, pp.68-71, 1994. 

[8] R. O. Grigoriev, M. C. Cross and H. G. Schuster, 
"Pinning control of spatiotemporal chaos," Phys. 
Rev. Lett., V 79, N 15, pp.2795-2798, 1997. 

[9] K.Zhou, J. C. Doyle and K. Glover, "Robust and 
optimal control," Prentice Hall, 1996, chapter 16. 

[10] G. E. Dullerud and S. G. Lall, "A new approach 
for analysis and synthesis of time varying systems," 
in Proc. 1997 IEEE/CDC. 



